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Abstract 

A nonlinear generalisation of Schro dinger's equation had previ- 
ously been obtained using information-theoretic arguments. The non- 
linearities in that equation were of a nonpolynomial form, equivalent to 
the occurence of higher-derivative nonlinear terms at all orders. Here 
we construct some exact solutions to that equation in 1 + 1 dimensions. 
On the half-line, the solutions resemble (exponentially damped) Bloch 
waves eventhough no external periodic potential is included. The so- 
lutions are nonperturbative as they do not reduce to solutions of the 
linear theory in the limit that the nonlinearity parameter vanishes. An 
intriguing feature of the solutions is their infinite degeneracy: for a 
given energy, there exists a very large arbitrariness in the normalisable 
wavefunctions. We also consider solutions to a q-deformed version of 
the nonlinear equation and discuss a natural discretisation implied by 
the nonpolynomiality. Finally, we contrast the properties of our so- 
lutions with other solutions of nonlinear Schrodinger equations in the 
literature and suggest some possible applications of our results in the 
domains of low-energy and high-energy physics. 
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1 Background 



Nonlinear extensions of Schro dinger's equation are often used as effective 
theories in domains such as optics [1] and condensed matter physics [2], 
the prototype being the Gross-Pitaevskii equation [3]. However, nonlinear 
Schrodinger equations have also been used to probe departures from exact 
quantum linearity jl] with several experiments setting upper bounds on some 
of the proposed nonlinearities [5] . The literature on theoretical investigations 
into nonlinear Schrodinger equations is immense, see for example the cita- 
tions in [6j. 

In Ref. [7] a particular nonlinear extension of Schrodinger's equation was 
motivated by information theory arguments [Bl E] similar to those used in 
statistical mechanics [ID] . For a single particle in one space dimension, the 
obtained equation was 

with 



F(p) = £ 



p(x + L) p(x) 



h 2 1 d'Jp 

H — (2) 

2m y/p dx 2 



and p(x) = ip*(x)ip(x). The energy scale £ is related to the length scale L 
through the relation 

< 3 > 

to ensure that the leading term in the small L expansion of (CQ) yields the 
usual linear Schrodinger equation. Note that if £ is finite, then L cannot be 
zero since h is non-vanishing. 

The nonlinear Schrodinger equation shares a number of properties 
with the linear Schrodinger equation. Of importance is that p(x) = still 
satisfies the usual continuity equation and thus has a sensible interpretation 
as probability. Furthermore, equation (CQ) still admits plane wave solutions 
ip = e tkx ~ lut and more generally, if ip is any solution of the equation, then so 
is \ip for an arbitrary constant A. 

As written, the nonlinear term F(p) is not invariant under the parity 
transformation x — > —x. A parity- invariant equation is easily obtained by 
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using a symmetrised form of Eq. (pQ) involving both postive and negative val- 
ues of L, 

F(p)^^(F + (p) + F-(p)) (4) 

where F + is the same as F but F~ uses — L instead of L. There is no problem 
in extending the nonlinear Schrodinger equation to more than one particle or 
higher dimensions, but in higher dimensions the equation is not rotationally 
invariant, leading to the suggestion [7J that spacetime symmetries might be 
linked to quantum linearity, with L possibly linked to gravitational effects 

m 

Given a nonlinear extension, it is natural to investigate whether any novel 
effects arise, such as the occurence of solitary waves, and to seek their experi- 
mental confirmation. In [7J some preliminary phenomenological consequences 
of the nonlinearity were suggested based mainly on qualitative and heuristic 
arguments. It was also noted in [7J that the equation had a unusual exact 
solution: Consider the free case,V = 0, of equation (jTJ) for a stationary state 
ip = e~ %Et l h (f)(x) with <f> real. Then substitution yields 



E 



h 2 d 2 



2m0 dx 2 
p(x) 



£ 



In 



p(x) 



£ 



In 



p(x + L) 
p(x — L) 



1 - 



p(x — L) 
p(x) 



+ 



h 2 d 2 ^p 
Im^fp dx 2 



(5) 



p(x + L) 



p(x) 



Clearly for any (f>(x) that is periodic in the nonlinearity scale L, 

(f)(x + L) = <j)(x), 



(6) 



(7) 



the right-hand side of the eigenvalue equation ([6]) vanishes and so E = 0. 
These zero-energy periodic solutions exist only for the exact equation ([6]) but 
not for perturbative approximations of the equation obtained by expanding 
to some power in L. For example, the perturbative approximation to (pO) 
using 



F(p) 



h 2 L 

4m 



+ 



p'p" 



+ 0(L 2 



(8) 



3p 3 2p 2 

where the prime denotes derivative with respect to x, does not support the 
class of periodic solutions of the exact equation mentioned above. 

Our purpose in this paper is to construct an even larger class of exact, 
localised, solutions for which the abovementioned zero-energy periodic solu- 
tions are a limiting case. We begin our construction in Sect. (2) for a single 



3 



particle confined to a box. A continuum spectrum of nonperturbative solu- 
tions is obtained which is not bounded below, even when a regularised version 
[7J of eq.([T]) is used. However when the equation is reconsidered on the half- 
line, then the spectrum of the solutions is bounded from below. In Sect. (3) 
we study a g-deformed version [7J of the equation ([I]) and again construct a 
continuum set of nonperturbative solutions whose spectrum is bounded from 
below for some values of the parameter q. 

Equation p]) suggests a obvious discretisation which we discuss in Sect. (4). 
Some approximate analytical solutions representing localised states on the in- 
finite line are also dicussed in Sect. (4): They might be useful seeds for further 
work. A discussion of potential applications and future directions appears in 
the final section. 



2 Regularisation and Solutions 

As the nonlinear term in ([1]) might have singularities where p(x) vanishes, it 
is convenient to consider a regularised version of that expression. Defining 

P±( x ) — p(x±r]L), (9) 

where the dimensionless parameter r\ takes values < r\ < 1, the regularised 
potential 



ln v , 1 0--v)p vp- 



7])p + 7]p + (1 — r])p + rjp + (1 — r])p_ + r\p 



(10) 



was constructed in [7j, and is to be used in ([2]) instead of the first term there. 
With this regularisation there are no singularities in the equation of motion. 
Formally, the rj parameter can be also used to interpolate between the fully 
nonlinear theory at r\ = 1 and the usual linear quantum mechanics at r\ = 0. 

We wish to find exact stationary states of the regularised nonlinear equa- 
tion. Notice that if 

p(x ± T]L) = g(L)p(x) (11) 

for some function g(L), then the term ffTUl) simplifies greatly, thus facilitating 
a solution of the nonlinear equation. Condition ffTTl) is equivalent to the 
eigenvalue equation for the translation operator, 

exp (±rjLd x )p(x) = g(L)p(x) , (12) 
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which is satisfied by a p(x) of exponential form. More generally, just as in the 
example for (jSJ), one may multiply that exponential solution by a function 
periodic in rjL and still obtain the simplification of ffTO]) . Thus the ansatzwe 
choose for our class of solutions is 

ip{x,t) =Cexp(— kx) a(x) exp (—iEt/h) (13) 

where C is the normalisation, k a real parameter and a(x) is any periodic 
function satisfying 

a(x ± r/L) = a(x) . (14) 

For the ansatz (1131) the regularised nonlinear Schrodinger equation re- 
duces to the following relation between the energy, E, and the other param- 
eters, 

(1 _^) = m + , (7 _ 1)1+ _i_ y (15) 

where 

7 = exp {—2KTjL) . (16) 
For definiteness one may take, 

a(x) = sin(^) , (17) 

then on the domain < x < Nt]L with N an integer, the normalisation 
factor is 

r -2 _ (1 ~ 1 N ) 16VT 2 

Ak (In 7 ) 2 + 16tt 2 ' 1 J 

With the choice ([17]) the wavefunction (flBl vanishes at the walls of the infinite 
well as required on physical grounds. Note that the equations (OQ [10]) require 
us to define the wavefunction simultaneously at points x, x ± rjL for each 
x in the domain [0, a]. This means that we must also give values to the 
wavefunction for some points outside the physical domain. The way we 
do this is to use ([TBI) throughout: This convention^] preserves not only the 
continuity of the wavefunction but ensures the validity of the eigenvalue 
equation (Tl5l) for all points in [0, a] . 



2 Of course one could adopt other prescriptions, such as periodic boundary conditions. 
If one imposed periodic boundary conditions on (|13p . then k = and one obtains the zero 
energy solution first found in Ref. [7]- That solution is still non-perturbative and highly 
degenerate, see (|T4|) . 
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Then, since on the finite x— domain the solution (I13II17P is normalisable 
for any — oo < k < oo, it means that for any given 77 the variable 7 can take 
values in the range < 7 < 00. Define 9 — 1 + 77(7 — 1). Then on the finite 
x— domain, (1 — 77) < 9 < 00 and so the function 

y = ln9 + ~ (19) 

is unbounded above, with a global minimum at 9 = 1 (corresponding to 
7 = 1). Thus the energy given by f|T5|) would be unbounded from below. 

Since it is necessary to have the energy bounded from below to ensure an 
energetically stable ground state, one clearly has to restrict the range of 9 in 
some natural way. This can be achieved by taking N — > 00 so that one is 
considering dynamics on the semi-infinite line < x < 00. Normalisability of 
the solutions now requires k > in (fT6|) . and so < 7 < 1 and (1— 77) < 9 < 1. 
Correspondingly, the energy is now restricted to 

-ind-,)-^). (20) 

As mentioned above, in the limit 77 — > the nonlinear equation reduces 
to the ususal linear equation. In the same limit, the energy ffTBT) becomes, 

ftV , . 

E — > , 21 

2m ' y ! 

where ([3]) has been used. However as fjl7H indicates, on the semi-infinite 
domain the wavefunction itself ( fl3l) does not reduce to that of the linear 
theory but instead undergoes infinitely fine oscillations. Thus even an in- 
finitesimal nonlinearity will leave its signature on the quantum ground state 
through nonperturbative effects: such possibilities are commonly studied in 
field theory contexts but we see them here at the quantum mechanical level. 

The above discussion has been for the parity- violating case of the nonlin- 
ear equation. If one considered the parity symmetric version of the equation 
using (j4j), then it is easy to check that even for dynamics on a semi- infinite 
region of the x-axis the energy would be unbounded below. 



3 A Q- deformation 

Another nonlinear Schrodinger equation that was constructed in [7] used a 
different information measure [12J, based on the deformed logarithm defined 
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by 



y 



9-1 



(22) 



q-1 

where the usual logarithm is recovered as the real parameter q — > 1. The net 
result is a nonlinear Schrodinger equation as in but with the first term 
in (121) replaced by the a q-deformed quantum potential 



Q 



(?) 

VL 
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q-1 




(23) 



with the other symbols having their previous meanings and q > in order 
to recover the usual linear Schrodinger equation in the small L limit. 

We can construct solutions to the q-deformed nonlinear Schrodinger equa- 
tion using the same ansatz (1X31 as before. The only thing that changes is the 
energy eigenvalue relation which now reads 



qE 

E q 



A 



(24) 



where we have defined 

A = exp (2«L) . (25) 

Again one finds that unless k is restricted in some natural way, the energy 
is unbounded from below. For dynamics on the half-line — oo < x < 0, 
normalisability of solutions requires k < and so < A < 1. Then one finds 
that the energy as given by (1241) is bounded from below (and above) if one 
chooses q > 1. 



4 Approximate Solutions and Discretisation 

The ansatz (Tl3j) does not permit localised solutions on the whole x-axis and 
it appears difficult to construct exact solutions with some other ansatz. How- 
ever, reasonably good approximate solutions can be constructed by patching 
the existing solutions as follows. Let 

27TT 

ip(x,t) = C'exp(-iEt/h) sin( ) (OixU+ix) + Oi-xUJx)) , (26) 

where 9{x) is the usual step function, C a normalisation factor and 

(f>±(x) = exp (tk±x) , (27) 



7 



with k± > 0. This normalisable and smooth wavefunction clearly satisfies the 
regularised nonlinear equation of Sect. (2) outside the small interval — 77L < 
x < rjL when only one patch of the wavefunction contributes. The values of 
n + and k_ are related to each other and to the energy E through a relation 
similar to eqs. fU5lll6l) . In the small interval around the origin eq. (1261) is not 
a solution so one requires a modification there. Thus the ansatz (12T)|) might 
be adopted as a trial for numerical searches of localised solutions, though we 
do not know if it would be useful in practice. 

In low-energy physics, nonlinear differential equations often arise as con- 
tinuum approximations to difference equations defining, for example, dy- 
namics on a lattice of atoms, see for example [13]. On the other hand, in 
high-energy physics, although spacetime is usually taken to be continuous, 
there have been numerous suggestions [H] that at a more fundamental level 
some form of discreteness should arise. Thus for both of those reasons, we 
look at discretisations of the nonlinear Schrodinger equations studied above. 
Let us focus on stationary states described by the basic equation §6§. It is 
remarkable that that equation has no explicit space derivatives, so the usual 
ambiguities of how to discretise a differential equation do not arise. In fact 
the structure of the nonpolynomial equation ([6]), which contains the function 
p(x) and its shifted values p(x±L), suggests a natural discretisation: Simply 
set L = 1, p(x) — > p n , p[x ± L) — > p n ±i, with n an integer, leading to a 
nonlinear difference equation for the energy eigenvalue. 

The above discretisation extends obviously to the regularised and q- 
deformed equations. The discrete forms of the equations might be more 
convenient for discovering new solutions in future investigations. 

5 Discussion 

The nonlinear Schrodinger equation ([I]), that was motivated by informa- 
tion theoretic arguments in [7], has a fundamentally different structure from 
other nonlinear extensions of Schrodinger's equation that we are aware of. 
Indeed, the nonlinear equation we have studied has the form of a differential- 
difference equation with points at x,x ± rjL simultaneously required in its 
definition. When expanded in L, the equation gives rise to higher derivatives 
of nonlinear terms at all orders. 

Furthermore, the nonlinearities in ([1]) are of a nonpolynomial form. It is 
interesting to note that nonpolynomial nonlinearities have been considered 
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in a physical context [12] but no higher order derivatives occured there. In 
a different context, nonlinear Schrodinger equations with derivatives at all 
orders have been considered in [13], but that nonlocality was only in the 
linear part of the equation. Thus the equation ([T]) combines features seen 
in separate equations in the literature, and perhaps it might find a physical 
application in those contexts. 

The nonlinear equation ([I]) has given rise to interesting solutions. Even 
in a bounded region of space, the energies of our class of solutions form a 
continuous spectrum and furthermore the wavefunctions are highly degener- 
ate, as indicated by the arbitrariness of the periodic function a(x). These 
features appear to be very different from various other exact solutions and 
discretisations of nonlinear Schrodinger equations studied in the literature, 
see for example [131 EE EEl EH] an d references therein. 

From a high-energy physics perspective, as originally envisaged in [7], 
the equation might be relevant for probing possible new physics at short 
distances. In this regard we note that higher dimensional, separable versions 
of the equation are available [7J. If one considers extra-dimensional scenarios, 
as is often the case in particle physics, the localised solutions on the half- 
line found here may be interpreted as extra-dimensional excitations that are 
localised near a brane by quantum nonlinearities. We also think that the 
result of Sect. (2) which indicated a nonperturbative quantum vacuum might 
have implications in cosmology. 

It is straighforward to construct some other exact solutions depending on 
more free parameters, for example by constructing the regularised version of 
the q-deformed equation in Sect. (3) or by replacing ffTT|) by an appropriate 
Fourier series. As for the equations we have considered, we have indicated 
some avenues for obtaining other localised solutions in Sect. (4). Finally we 
note that the solutions studied in this paper do not deform continuously to 
solutions of the linear theory as the nonlinearity vanishes. The perturbative 
effect of the nonlinearity on the usual solutions of the linear Schrodinger 
equation are studied elsewhere [IS] . 
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